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Abstract There is a growing population of relativistically relevant minor 
bodies in the Solar System and a growing population of massive extrasolar 
planets with orbits very close to the central star where relativistic effects 
should have some signature. Our purpose is to review how general relativity 
affects the orbital dynamics of the planetary systems and to define a suitable 
relativistic correction for Solar System orbital studies when only point masses 
are considered. Using relativistic formulae for the N body problem suited for 
a planetary system given in the literature we present a series of numerical 
orbital integrations designed to test the relevance of the effects due to the 
general theory of relativity in the case of our Solar System. Comparison 
between different algorithms for accounting for the relativistic corrections are 
performed. Relativistic effects generated by the Sun or by the central star 
are the most relevant ones and produce evident modifications in the secular 
dynamics of the inner Solar System. The Kozai mechanism, for example, is 
modified due to the relativistic effects on the argument of the perihelion. 
Relativistic effects generated by planets instead are of very low relevance but 
detectable in numerical simulations. 

Keywords celestial mechanics • relativistic effects • Solar System • small 
bodies 



1 Introduction 

General Relativity Theory (GRT), while widely used in some areas of astro- 
physics, has not yet been fully exploited in the study of planetary system 
dynamics, despite the fact that one of the most famous relativistic effects, 
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that of the precession of perihelia, shows that Relativity can play a role, 
albeit somewhat secondary, in the evolution of planetary systems. 

The design of a detailed general relativistic reference frame for use in the 
description of our Solar System was achieved in a compelling though very 
technical way b y Brumberg and Kopeikin and also by Damour et a l. in a series 
of publications (iKopeikin Il988l: iBrumberg fc Kopeikinl . Il989al lbl: iBrumberel 
Il99lt iDamour et all ll99lL Il992l . Il993h . In those papers, each object in the 



system was equipped with a local reference frame (for measuring internal de- 
grees of freedom) which was fitted within a global reference fram e. This work 
culm inated in resolutions of the IAU 2000 General Assembly (jSoffel et all 
2003). The methods of Brumberg and Kopeikin, and Damour et al., allow 
for the development of relativistic equations of motion for extended objects, 
but for purposes of studying the general dynamics of planetary systems, the 
limit of point (monopole) masses can be expected to give good qualitative 
results. 

Also, there has been exhaustive studies of the validity of GR T by means 
of the so called Parameterized Post Newtonian ( PPN) theories ( Willi . Il98ll 
119931: iKlioner fc Soffell . feOQOl ; iKopeikin fe Vlasovl . 12004 ) and co mpari s on wi th 
observations made in the Solar System (for a recent review, see IWilll (2006)). 
Nevertheless, the effects measured for that purpose are in general not related 
to any dynamical study of the planetary system. There have also been a 
number of studies of Post Newtonian effects, such a s the Lens-Thirr ing effect 
i n the motion of satellites (see for example lloriol ( 20051 ) and also ICiufolinil 
(2007)), but again these were not concerned with planetary dynamics studies. 

Except for very precise simulations (for example iQuinn et al. (1991). 
IVaradi et al.1 (|2003l) l relativistic effects are in general not taken into account 
in the study of Solar System dynamics, and when they are, only the effects 
due to the Sun are considered, using a simplifica tion of the above mentioned 
more complete PPN model found in I Willi ( 198ll ). This simplified approxima- 
tion was also used in the study of some extrasolar planetary systems. For 
example. iBarnes fc Quinnl ( 200ll ) studied the s tability of the i>-Andro medae 
planetary system using that approximation and lNagasawa fc Linl(|200Rh showed 
strong differences between classic and relativistic evolutions of planet b of 
that system. A more general study of the contribution of t he relativistic terms 
to the secular evolution of extrasolar systems was done by[Ad ams fc LaughlirJ 
(|2006h . In the case of extrasolar planets orbiting Pulsars, a specific s t udy o f 
relativistic dynamics is still pending, as stated in (Gozdzwies ki et aill2005l ). 

In the case of minor bodies, as explained by Shahid- Saless &: Yeomansl 
( 1994 ) GRT should be used when studying low perihelia orbits. They also 
pointed out that inconsistences arise when using newtonian dynamics with 
masses derived from relativistically derived ephemerides. 

Apart from the aforementioned studies, relativistic effects were included 
in the pioneering work at high precision ephemerides at the JPL, and by now 
they have a measurable impact o n observ a tions, with today's high precision 
astronomical techniques (see e.g. iPitieval (|2005l ) for recent advances in the 
numerical integrations of planetary ephemerides). 

It seems without doubt that future development of dynamical astronomy 
and experimental gravitation will fundamentally rely on the study of General 
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Relativity effects. This is a comm on point of view ex pressed in resolutions of 
the IAU 2000 General Assembly (jSoffel et al.1 . l2003h 

The aim of this paper is to evaluate how GRT affects the orbital dynam- 
ics of the Solar System, in the limit of point masses, by means of a series 
of numerical tests. By extension we can extract some conclusions for the ex- 
trasolar systems. Hopefully, the numerical studies here presented will help 
to elucidate in which circumstances it is necessary to account for relativistic 
effects. 



2 The Post-Newtonian Algorithm 

We are interested in the first corrections to the classical Newtonian equations 
of motion of a system of point masses, as derived from the GRT. They can 
be found by means of an approximation of this last theory, known as the 
Post Newtonian approximation. This approximation can be accomplished by 
making an expansion of the GRT equations in terms of v/c, where the v 
are the velocities, and c is the speed of light. Equation of motion can thus 
be found for point masses, which include the Newtonian term, as well as 
a number of other terms which are suppressed by factors of 1/c. T he PN 
equation of motion used in this work is given bv lNewhall et aL 
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with fj,i standing for Grrii, Vj is the barycentric velocity of rrij and dots 
denoting usual derivation with respect to an inertial time, making ?-,- the 
acceleration of rrij . Here, "velocity", and "acceleration" must be understood 
in a coordinate sense. These quantities do not in fact correspond exactly to 
the velocity and acceleration as measured by any physical observer - a subtle 
difference not found in Newtonian Mechanics. The summation is over all mas- 
sive bodies including the Sun and the origin is the barycenter of the system. 
This expression is known as the Einstein-Infeld-H offman equation, and the 
firs t version of it is fo und in lEinstein et al.l (|l938l ) . For a modern deduction, 
see Brumberg ( 2007 ). The expression also coincides with the one given by 



iDamour et al 



19911 ) in the case of point monopole masses. A pedagogical 



derivation of this equation, starting from GRT, is given in appendix [Al 

As it stands, this algorithm would greatly slow down any numerical calcu- 
lation, and it is only useful in situations where the relativistic contribution of 
every object in the system is to be taken into account. For the simplified case 
where only one object (e.g. the Sun) contributes with relevant corrections, 
the sums in {1} are changed into a single term, and in the reference frame of 
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this central, massive object with mass M, the relativistic correction would 
read 



A .. GM 

Ay = -^r 



AGM 2 , 
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(2) 



which is the correcti on proposed by I Anderson et al.1 ( 1975( ) (in fact, the 



equation as written in lAnderson et al.l ( 1975 ) depends on two PPN pa- 
rameters (3 and 7, which are equal to unity in the G RT case) an d used 
by most papers that consider relativistic correcti ons (|Quinn et all 1 19911: 
IShahid-Saless fe Yeomaiil . ll994lVaradi et al.l . l2003l) . This algorithm, although 
less exact than Eq. fTJ), is computationally much more affordable and, as we 
will show, allows for the most important effects. Anyway, in certain cases the 
complete algorithm fl} may become necessary, when there are more than 
just one relativistically relevant object in the system under consideration. 
Sometimes, in very simplified analysis, the second term in the right hand 
side of (|2|) is ignored but this can only be justified when considering near 
circular orbits. 

In going from equations (JTJ) to (J2J) one is not taking into account possible 
relativistic effects due to the changes in the velocity of the central object 
caused by other, large but not necessarily relativistically relevant, objects. 
In the case of the Solar System, the influence of Jupiter on the motion of 
the Sun could in principle introduce such an effect, but a simple order of 
magnitude analysis shows that this effect is highly suppressed and can be 
consistently taken out. 

It is possible to show that the only secular contribution -that means non 
short period- from relativistic effect generated by the central star affects only 
the argument of the perihelion uj and the mean anomaly M which is related 
to an effect on the time of perihelion passage as measured from a far observer 
at rest. The effect on u> due to the relativistic effects of the Sun is giv en by 
Auj = 0.0384/ (o 5 / 2 ( 1 - e 2 )) arcseconds per year whe re a is in AUs (jWilll . 
119811 : ISitarskil 119831 : IShahid-Saless fc Yeomani . Il994h . For low eccentricity 



orbits the variation in M as measured by a far observer at rest is given 
approximately by AM ~ — 0.115/(a 5 / 2 v / l — e 2 ) arcseconds pe r year and i t 
is possible to show that the short period effects are neglig ible (lIorioLl2005l) . 



Both secular effects should be considered, for example, in the case of comets 
when obtaining their non gravitat ional forces from the observed temporal 
evolution of the orbital elements ( Yeomans et a l.. 2005). These analytical 



predictions due to pure relativistic effects can be easily verified by means 
of a numerical integration of a masslcss particle orbiting the Sun under the 
effects of Eq. ([2]). We show as an example a numerical integration of a particle 
with orbital elements similar to 2P/Encke at Fig. [l] The short period effects 
are undetectable if any at all, and the effects on uj are exactly as predicted. 
The effects on mean anomaly M are not exactly as predicted because of the 
high value of e. The other orbital elements remain unchanged. In the figure 
the time is as measured from a distant observer at rest. 

Some symplctic numerical integrators use an artificial separable H amilto- 

nian that can mimic the effects on uj and M that Eq. ((2]) generates ( Saha fe Tremairiel 



119921 ). Using this artifact the relativistic effects can be included in the per- 
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turbing Hamiltonian and sympletic form is maintained but this approxima- 
tion, usually used in studies of our planetary system, is only valid for nearly 
constant and not very high eccentricity orbits. 

We are interested in exploring how relativistic effects generated by Eqs. 
([T]) and {21 affect the dynamical evolution of our planetary system and the 
dynamics of small bodies. In order to analyze the relativistic effects in the 
Solar System we implemented three different models: one using classic new- 
tonian accelerations (model C), another which considers only the relativistic 
effects generated by the Sun (Rs) and the last with the expression that takes 
into account the relativistic effects generated by all massive bodies (R a «). 



3 Effects in Our Planetary System 

Our planetary system is mainly composed by an inner set of terrestrial low 
mass planets and an outer set of giant very massive planets. As they evolve in 
different timescales we will analyze both two groups separately. Our model 
for planetary system includes the Sun and the eight mayor planets from 
Mercury to Neptune and we consider the Earth-Moon system as unique body 
lo cated at its baricen ter. This model is simpler than, for example, the one 
of iQuinn et all (|l99ll) so our results will not be more precise. Our aim is to 



study the relevance of the relativistic corrections in order to have an idea if 
they are or not important in modeling the orbital dynamics of bodies evolving 
in the Solar System. 



3.1 The inner planetary system 

The inner planetary system is strongly perturbed by the giant planets and 
it is known that the relativisti c terms due to the Sun have sign atures in 
timescales of millions of years (lLaskarl Il988l [Q mnn et all 1 199 111 . Are the 



relativistic effects generated by the planets (model R a /;) relevant or on the 
contrary are the Sun's relativistic effects (model Rs) enough to explain the 
dynamics of the inner planetary system? 

In order to investigate this point we have integrated all the planetary 
system from Mercury to Neptune by 2 millon years. We performed three 
different runs accordin g to the three di fferent models: 1) model C using 
MERCURY integrator (|Chambersl 1 19991) with the Bulirsch-Stoer method 



including only classical newtonian gravitation, 2) model Rs using the same 
integrator with BS method including an user defined acceleration given by 
Sun's relativistic effects according to Eq. (|2|) and 3) model Knu wit h our 
own fully relativistic routine implemented using RA15 (|Everhartlll985h that 
includes relativistic effects from all massive bodies according to Eq. (JTJ) . 
Runs with models C and Rs were repeated for control using the classi c 
and relativistic versions of the integrator EVORB ([Fernandez et all (|2002l h 
|http:/ /ww w .fisica.edu.uy/~gallar do /evorb.html| ) . 

Figures show the evolutions of the pla nets from Me rcury to Mars. 
It is possible to see that the main results from lLaskarl (|19881 ) are recovered. 



Notable differences between classic and relativistic models are evident in the 
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case of Mercury (Fig. [2]). In the case of Venus (Fig. [3j) and the Earth (Fig. HJ, 
the relativistic integrations present also some variations with respect to the 
classical one in the evolution of eccentricities and inclinations, in addition to 
a precession of the perihelion. For Mars, effects are of much lower importance 
and start to be appreciable at the end of the integration, also mainly in the 
evolution of eccentricity and inclination. 

For the considered time intervals the integration results from both rel- 
ativistic models R5 and R Q ;/ were basically undistinguishable. We want to 
stress that both relativistic models Rs and R a /; were implemented in very 
different numerical integrators with very different algorithms (BS and RA15 
respectively), so the coincidence of the results is a prove of their precision 
and a confirmation that the differences detected with model C are real and 
not a numerical artifact. 

We tested the relevance of the relativistic effects by running an integra- 
tion that included Earth plus the outer planetary system (that is without the 
less massive terrestrial planets) for a time span of 2 millon years. By doing 
this, we eliminate the forced modes generated by the other terrestrial plan- 
ets and we iso late the relativistic effects from the inner Solar System chaos 
(Laskar, 1996). The resulting evolution of Earth's eccentricity and specially 
the inclination are much less affected by the relativistic corrections. Then, 
the discrepancies between classic and relativistic evolutions found in the real 
Earth in Fig. |4] are mainly due to the complexity that characterize the dy- 
namics of the complete inner planetary system, making the small relativistic 
effects on the inner planets to interact in an complex way. 



3.2 The outer planetary system 

We have shown that relativistic effects generated by the Sun are relevant but 
the ones generated by the planets are negligible in the inner planetary system. 
Can we conclude the same for the outer more massive planetary system? 
We are interested in relativistic effects due to the massive outer planetary 
system and to compare with the relativistic effects due exclusively to the 
Sun, that means model R a ;/ versus Rg. In this experiment we integrated 
the outer planetary system from Jupiter to Neptune by 401 millon years. We 
performed three different runs according to the three different models defined 
earlier. 

In Figs. [6] and [7] we show the last years of the evolution of the eccentricity 
and inclination respectively for the four giant planets. There is a small shift 
in e and a very small shift in i between the di fferent model s which should 
affect very slightly the fundamental frequencies ( Laskarlll988f ) but the global 
behavior is not substantially changed. 

Results for the more massive and near Jupiter and Saturn are somehow 
different than results for the less massive and remote planets Uranus and 
Neptune. Analyzing the eccentricity showed at Fig. [6] we can see that model 
R a z; departs from R5 and both from the classic model C. In Uranus both 
relativistic models are almost coincident and they depart from C whereas 
for Neptune differences are something small. Analyzing the inclinations for 
Jupiter and Saturn in Fig. [7] we can see that model C and Rs coincide and 



7 



they depart from R a /;, so mutual relativistic effects between Jupiter and 
Saturn seem to be more important than the relativistic effects of the Sun on 
Jupiter and Saturn. No differences in inclination are appreciated in Uranus 
and Neptune. We can interpret this as follows: the two most massive and 
near planets Jupiter and Saturn generate some additional perturbations in 
model R Q ;/ but these effects are not enough in the more distant and less 
massive planets Uranus and Neptune. Effects in the inclination are of lower 
magnitude than in the eccentricity and only in the more massive and near 
planets R a /; show up. 

It is known that quasi resonant angles like 2Aj — 5As or 2Ajy — \u 
where A is the mean longitude can generate non negligible perturbations on 
minor bodies (iFerraz-Mello et all . Il998t iMichtchenko fe Ferraz-Mellol 1200 it 
M arzari fc Scholll . 12002t h A modification of the time evolution of the quasi 
resonant angles due to relativistic effects could introduce dynamical effects 
in the time evolution of these minor bodies. Analyzing our integrations we 
conclude that the circulation periods of both critical angles remain practi- 
cally unchanged by the relativistic terms, so that there do not exist relevant 
relativistic effects on these angles. 



4 Effects in Solar System's Low (a, q) Orbits 

4.1 Secular dynamics of fictitious particles 

According to the expression for the relativistic effect Acu we can infer that 
the low (a, q) orbits should be the most sensible to the relativistic effects, and 
according to the results of section 3.1, if we exclude planet crossing orbits, 
we can ignore the relativistic effects generated by the planets in the case of 
low (a, q) orbits. So, we have integrated by 10 5 years the system composed 
by the Sun, the four giant planets and a population of test particles with low 
perihelion (0.05 < q < 0.5 All) and low semimajor axis (0.2 < a < 2 All) 
with and without the relativistic terms due to the Sun, that is, models Rs 
and C respectively. We also performed a run with the model R5 including the 
terms corresponding to the solar quadrupole moment of the Sun J2 = 2x 10~ 7 
(|Pireaux fc Rozelotl I2003T) obtaining a completely negligible effect, so these 
effects can be ignored in comparison with the relativistic effects due to the 
Sun. 

We considered the time evolution of the elements e,i,uj and computed 
the differences between both runs, classic and relativistic. The maximum 
differences Aoj, Ae and Ai found during the integration are showed at Fig. 
[5] as a function of the initial semimajor axis and perihelion distance. That 
figure was constructed with 100 test particles all with initial i = 20°, J? = 
0° , uj = 60° , so it is not a map valid for any set of orbital elements but give 
us an idea of how relativity modify the secular dynamics in that region. 

This experiment designed excluding the terrestrial planets avoids the 
chaotic dynamics typical of objects with close encounters with them and 
the differences between both runs give us and indication of the differences 
generated by the relativistic terms in the secular, more regular, dynamics. 
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B esides some secular resonan ces present in this region ([Michel fc Thomad . 

119961 : iMichel fc Froeschlel Il997h the secular evolution that dominate this re- 
gion is the Kozai mechanism which produces strong variations in (e, i) linked 
to the evolution of ui maintaining a constant value of H = \/l — e 2 cos i, 
where the inclination is measured with respect to the invariable plane of the 
Solar System. The relativistic terms affect lo mod ifying the Kozai mechanism 
and in consequence the time evolution of (e,i) (|Morbidellil . l2002t ) . 

At Fig. [S]we show an example of how relativity modifies the Kozai mech- 
anism. It is a typical secular evolution with constant a — 0.4 AU in both 
models within all the time integration period. Strong differences between the 
other orbital elements appear conserving H nearly constant. 

In the real Solar System these effects can be increased several times due 
to the secular forced modes generated by the terrestrial planets but in the 
case of orbits having close encounters with the planets the chaotic dynamics 
will overcome any improvement to the model introduced by the relativistic 
corrections. 



4.2 The known relativistic asteroid population 

After the successfully explanation of the precession of perihelion of Mercury 
by Einstein's GRT, an unavoidable step was to st udy the evolut i on of the first 
very low perihelion asteroid known, 1566 Icarus ( Shapiro et al.lfl968h . where 



relativistic effects were expected and observed in its dynamical evolution. 
Icarus, being maybe the most studied "relativistic" asteroid, is by no means 
the only one to deserve such treatment. Mercury is the Solar System's body 
with greatest relativistic Au> (0.4304 "/yr) but there are at present around 
half hundred asteroids with Alo greater than Icarus (Au = 0.1006 "/yr ). The 
objects with the greatest values for Acu are shown at Table Q] which was con- 
structed from the Near Earth Asteroids database of the Minor Planet Center 
and ca n be considered as an update of Table 1 of IShahid-Saless fe Yeomansl 
(1994). These asteroids are the ones with the most important relativistic ef- 



fects on uj but not necessarily the ones with the most diverging evolutions 
from classic model C. Divergence depends also on the degree of chaos of their 
dynamics. 

We have integrated the full planetary system plus all these asteroids 
and also some others for a time span of 10 5 years using both the classi- 
cal and relativistic models C and R5. All asteroids from the table exhibit 
chaotic evolutions with several close encounters with terrestrial planets that 
generate quick d epartures between classic and relativistic trajectori es. Res- 
onance sticking ( Duncan fe Levisonl . 1 1997t iLvkawka fc Mukal . l2006h is also 



pre sent with s witch es between resonances that can be easily identified follow- 
ing [Galardc] (|2006[ ). Asteroid 2004 XY60 is a good example (Fig. [TO]) of an 



evolution dominated by close encounters plus resonance sticking. It starts in- 
side resonance 6:5 with Venus and after 13000 years the resonance's strength 
drops and high order resonances with Earth start to dominate following dif- 
ferent dynamical trajectories the classic and the relativistic models. 

As an extreme case we can mention 2000 LK, not included in Table [TJ in 
which the relativistic particle hits the Sun after 83000 years, and the classical 
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one survive during the whole integration (Fig. [TT|). A more interesting case 
is 2003 CP20 (Fig. [T2|) neither included in the table, which did not have any 
close encounter with planets in any of the runs but it is strongl y affected 
by ext erior resonance 27:28 with Venus as we have found following lGallard"ol 
(l2006h . The critical angle a — —27\y en + 28A — m is shown at Fig. [TBI for 
the two models C and Rg. In this case the Kozai mechanism is also present 
generating important variatio ns in e, i in analogy t o other situations observed 
in the transneptunian region ([Gomes et all l2005h . 

We have also studied some real asteroids and fictitious particles in mean 
motion resonances with Jupiter obtaining very different evolutions between 
models C and Rs. This is an expected result because the chaotic nature of 
the resonant dynamics makes that small variation in the models to generate 
exponential divergences. For completeness we have studied the dynamical 
evolution of fictitious particles in the transneptunian region but no relativistic 
effects were found there. 



5 Extrasolar Systems 

The relevance o f relativistic effec t s in e xtras olar systems were alread y shown 
by, for example, Mardling fc Linl (|2004 ) and lNagasawa fc Lid (|2005h . In par- 
ticular iNagasawa fc Linl (|2005f) showed that relativistic t erms have dramati- 
cally a ffected the evolution of planet b in v Andromedae. lAdams fe Laughlir] 



( 20061 ) have extended the secular interaction theory including GR terms and 



they found that relativistic corrections are relevant for systems with close 
planets (~4 days orbits). Then, it is not necessary to insist here on the im- 
portance of including relativistic corrections specially when semimajor axes 
are of the order of 10 _1 AU or less, but we can analyze the relevance of the 
use of model R a u instead of Rs • 

With this in mind we integrated a fictitious planetary system with low 
eccentricity and close orbits of massive planets (see Table [5} using the three 
models. We avoided experiments with very close orbits in order to avoid 
chaotic dynamics. The global dynamics obtained by the three models are very 
similar and regular due to the small eccentricity of the orbits. Differences are 
mainly in the frequencies of the secular evolution. We have fourier analyzed 
one of the fundamental frequencies of the system obtaining / = 1.63677 x 
10~ 3 yr- 1 for model C, / = 1.63682 x 10~ 3 yr" 1 for model R s and / = 
1.63681 x 10~ 3 yr _1 for model R a /;. We have checked that the differences were 
not a numerical artifact repeating the simulation R a ;/ with lower precision 
parameter for RA15 and verifying that the new frequency was exactly the 
same. Then, differences between models Rs and R a ;z seem to be real and 
not due to numerical artifacts but, at least in this example, negligibly small. 
We repeated the experiment with the same planetary system but in this case 
with smaller semimajor axes (0.2 and 0.8 AU) and an analogue behavior was 
observed regarding the differences between fundamental frequencies in the 
three models. 

However, relativistic effects due to the planets (model K a ii) could be 
important for planetary systems in mean motion resonances or under the 
mechanism of perihelion alignment or ant i- alignment evolving very near the 
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separatrix in the phase space because a small difference in the evolution of w 
could generate a very different behavior in the phase space with an associated 
very different evolution in e. 

The problem of the relevance of model R a ii deserves a more detailed study. 
In principle we can conclude that in planetary systems, even with low e, it 
is worthwhile to consider rclativistic corrections at least with model Rg and 
a test with model R a ii will be convenient. 



6 Conclusions 

The secular relativistic effects generated by the Sun are relevant for the argu- 
ment of the perihelion and mean anomaly (related to the time of perihelion 
passage) whereas the short period effects are vanishingly small. In our plane- 
tary system relativistic effects generated by the Sun are appreciable specially 
in the inner Solar System, but those generated by the planets are very small 
and are only detectable in the long time scales in the giant planets, specially 
Jupiter and Saturn. 

Due to the secular dynamics the departure between classic and relativistic 
trajectories is enhanced in the inner Solar System. In particular the Kozai 
mechanism is substantially modified driving notorious differences in the evo- 
lution of the eccentricity and inclination. 

Nowadays there exist a noticeable population of rclativistic minor bodies 
in the Solar System. Around a half hundred asteroids at present exhibit 
greater relativistic effects than 1566 Icarus. The relativistic corrections due to 
the Sun are enough for the correct modeling of the dynamics of low perihelion 
or low semimajor axis bodies. Minor bodies evolving having close approaches 
to the planets or being in resonances are the most affected by the relativistic 
corrections not because of the magnitude of the relativistic corrections but 
due to the chaotic nature of their dynamics. Relativistic corrections should 
also be applied to comets but in this case in order to be consistent a suitable 
model for the non gravitational forces is necessary. 

In general, in extrasolar planets the relativistic corrections due to the star 
are enough and, in some cases, are indispensable as in the case of systems 
evolving near a separatrix of trajectories in the phase space. In some partic- 
ular cases the relativistic effects generated by the planets could have some 
importance, but this point deserves future research. 

In the case of exoplanets orbiting around pulsars, a more detailed study 
needs to be performed. Such a study could be done following the same lines 
as those presented in this paper, and it constitutes the obvious direction for 
future research. As for the case of binary pulsar systems, its study is more 
complex and the a lgorithm s used in this paper would need to be improved. 

Some years ago lSitarskil ( 1983T ) concluded "...it seems that in all the mod- 



ern investigations it is the very time to replace Newtonian equations of mo- 
tions by those resulting from general relativity theory" . Following the results 
here presented, it seems that relativistic corrections are not unavoidable for 
all dynamical studies, but they are necessary for the precise dynamical mod- 
eling of bodies evolving in the inner Solar System and for the dynamical 
modeling of close extrasolar planets. 
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A Appendix: A Closer Look at the Post-Newtonian Algorithm 

General Relativity is a theory written in a language quite unlike the one of newto- 
nian gravitation. In GR, space-time is equipped with a metric <? M „ , which measures 
distances and angles at every point. There does not exist any gravitational force 
in the theory, and the gravitational interactions are given by the form of the met- 
ric itself. Free particles in the theory follow geodesies, that is, trajectories which 
minimize distances as given by the metric. 

Nevertheless, as any new theory, GR must, and of course does, recover the 
classical theory as a limit of weak gravitational fields and low velocities. This is 
what is called the newtoman limit of the theory, and in it the classical field <j> is 
recove r ed as a term in the te mporal part of the metric g^ v of space— time IjWeinberd 
(p7l ); IThorne et alJ <fl973 )) 

300 = -(1 + 20) 

In our case, we take this limit and improve it using an expansion in the typical 
velocities of the particles of the system v 2 ~ GM/r , following IWeinberel (|1973 ). 
For velocities which are small compared with the speed of light c, that means 
GM/rc 2 <C 1, this approximation is sufficient to describe the dynamics of point 
masses in the Solar System. Following the usual convention in GR, we work in units 
where c = 1. We can then write 



goo = + +Soo 3 +••• 
9H = Sij + glf + gif + ... 

(3) , (5) , 

9io =9h +9h +■■■ 

where numbers in parenthesis denote the order in v/c of the terms, and the Latin 
indexes represent the spatial components. Notice that, in order to preserve prop- 
erties of symmetry under time reversal at this order of the expansion (which is a 
valid symmetry, due to the fact that gravitational radiation e ffects are of a higher 
order in the expansion (IWeinberel [1973 : ITS orne et all 119731 ). only terms with a 
given parity are allowed in the expansion. 

Similar expansions can be produced for the inverse metric tensor <7 M " and for the 
Ricci tensor of curvature R a /3- In GR, curvature of space-time is given by second 
derivatives of the metric tensor, and measures how much geodesies are different 
to straight lines. The relation between these expansions can be readily deduced, 
giving expressions in the form of 

{2) _ 1 1 9^ _ 1 _ 1 % w 

iJ 2 dx*dxi 2 dx^xi 2 dx k dx* 2 dx k dx i ^ 

for the spatial components, making use of Einstein's summation convention. 

These expressions can be simplified by the use of some arbitrary gauge condi- 
tions, which are permitted due to the invariance properties of GR theory under 
arbitrary coordinate transformations. This kind of freedom i s analogou s to th ose 
that exist in electromagnetic theory. In our case we followed Weinberg (|1972f ) in 
choosing the so-called Harmonic Coordinate Conditions, which force the metric to 
satisfy 



r x = sTl^, = 
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where r£ v are the Christoffel symbols, related to the first derivatives of the metric, 
and also to the parallel transport of vectors along a curve in space-time. The use 
of these conditions impose restrictions to our expansion terms, e.g. 



0: 



1 dgfp dg. 

2 dt 



_1\ IM 

dx i 2 dt 



(2) 



these restriction help to simplify our equations. 

Then we can proceed to make use of Einstein's equation, which relates in each 
point the curvature of space-time to the presence of matter and energy: 



7^g^T\) 



with 



where T M „ is the energy-momentum tensor, including the density of energy in space- 
time. In the case of point particles, this density will reduce to a summation over 
Dirac's 8 functions at the position of each particle. This tensor can also be expanded 
in terms of v/c, and when we substitute the expansions in Einstein's equation we 
reach equations of the form 



v 2 ^ 



V 2 g$ = IGnGT 



_ 87rGT (o)oo 

+2 ' "ij 



Of® 



dt 2 aiJ dxtdxi \dx* \ dx* 
SttG [t< 2 > 00 - 2g$TW°° + T< a >«] 



(l)iO 



-8nG6i d T 



(0)00 



And we find, as was to be expected 
g$ = -20 = 2G 



T 



(0)00/ / 



where T (0)0 ° is just the mass density of the system, and <f> the classical potential. 

We can now define some new potentials, which will play the role of corrections 
to the classical theory. 



o (3) 
9io 







where 



Ci(r,t) = -4G [ r (1)l0 ^4-dV 
J r - r 



g$ = -20 2 - 2</> 



where 



V(r,t) 



/ 



2 Ur',t) 



+ GT (2)00 (r',t) + GT {2) "(r',t) 



4tt dt 2 



d r 



It is possible now to write a relativistic equation of motion in terms of these 
potentials. As we have already stated, in GR test particles move following geodesies, 
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and the equation for a geodesic is given, in terms of a generic curve length parameter 
T by 

dr 2 vX dr dr 

where the Christoffel symbols, which are needed for making covariant derivatives in 
a curved space-time, are also written in terms of the metric, and of its expansion. 
Splitting up of temporal and spatial indexes leads to the expression 



dt 2 



7~i2 r» j-ii dx j dx dx 

-r 00 -2r oj — -r jk — — 



r u + 2r; 



dx>_ 
It 



+ r 



dx j dx k 
~dT~dT 



dx i 
dt 



finally, after recovering the factors of c for a more intuitive expression, substitution 
of the three potentials give 



Given this equation of motion, it becomes possible to understand the physical 
meaning of the three potentials, cj>, ip and The case of (j> is easiest, as it is nothing 
but the classical gravitational potential from Newton's theory, thus recovering the 
classical limit when v < c. The vector potential ( is quite similar to the vector 
potential A in electrodynamics. Just as the magnetic field can be seen as a (special) 
relativistic effect of the electric field, so £ appears as a relativistic correction to the 
classical field, which exerts a velocity-dependent force normal to the distance vector 
and to the relative velocity of two massive objects. This field is sometimes known 
in the literature as the gravitomagnetic field, for its obvious analogies with the 
magnetic field in electrodynamics. 

The ip scalar field contains corrections to <f> due to many relativistic effects. 
It can be shown that the time derivative term in the ip defining integral accounts 
for the delay of the gravitational signal, which travels at the speed of light, and 
which moves following a geodesic instead of a straight line. The other two terms 
are corrections due to the equivalence between mass and energy, by means of which 
the kinetic and potential energies of the particles modify the classical field. 

The equation of motion contains the classical equation of Newton as its first 
term, and many small corrections to it. The equation in itself can prove to be very 
insightful for a theorist, but it is not fit yet for use in a numerical algorithm. We 
would like to find a more straightforward algorithm, in which the acceleration of 
each point particle would be given explicitly by the positions and velocities of every 
particle in the system. 

To accomplish this we must use the relativistic expression for the momentum- 
energy tensor of a system of particles, in order to find <f>, %p and £ as functions of 
positions and velocities. In the case of ip, the delay-of-signal term can prove very 
hard to simplify. Some methods and approximations from electrodynamics can be 
used, as the d elay-of-s i gnal t erm is almost identical to the Lenard-Jones electric 
potential (e.g [Griffiths ( 1999)). After some algebra it is possible to arrive to the 
formula given bv lNewhall et al.l |l983), which is more explicit - and numerically 
useful - albeit somewhat less physically meaningful: 
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, v ? 4 „ ^ 3 /(^-r.O-vA 2 , 1 



+ ^ E ^ [(" - ■ - 3vi)l (v* - v 3 ) + ^ ]T ^ 
which is the same as eq. (JTJ) in the main text. 
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Table 1 Solar System's bodies most affected by relativistic effects generated by 
the Sun. 



Designation 


a (AU) 




e 


V / 


Auj 


C'/yr) 


Mercury 


0.387 





206 


7.0 





430 


2004 XY60 


0.640 





797 


23.7 





321 


2000 BD19 


0.876 





895 


25.6 





268 


2006 CJ 


0.676 





755 


10.2 





238 


(66391) 1999 KW4 


0.642 





688 


38.8 





221 


1995 CR 


0.907 





869 


4.0 





200 


1999 MN 


0.674 





665 


2.0 





185 


2006 KZ39 


0.616 





525 


9.3 





178 


2001 TD45 


0.797 





777 


25.4 





171 


2004 JG6 


0.635 





531 


18.9 





166 


1998 SO 


0.731 





699 


30.3 





164 


(85953) 1999 FK21 


0.739 





703 


12.5 





162 


2004 XZ130 


0.618 





455 


2.9 





161 


2006 JF42 


0.672 





582 


5.9 





157 


2005 WS3 


0.672 





576 


23.0 





155 


2004 UL 


1.266 





927 


23.7 





151 



Table 2 A fictitious planetary system adopted to evaluate the effects of models 
C, Rs and K a ii- The mass of the central star was taken equal to Mq. 



Object 


a(AU) 


e 


i( u ) 


m/m.jup 


Planet b 


0.6 


0.04 


1.04 


1.0 


Planet c 


1.0 


0.02 


1.45 


1.0 



0.06 




-0.02 1 1 ' ' 1 1 

2000 4000 6000 8000 10000 

time (years) 

Fig. 1 A numerical integration considering only the Sun and a massless particle 
with orbital elements similar to 2P/Encke (o = 2.22, e = 0.85) designed to show 
the relativistic effects on lu and M due to Eq. ([2| • Variations are Auj 
and AM = M re i — M c \ a . The short period effects are undetectable. 
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Fig. 5 Mars' evolution. Top: eccentricity; bottom: inclination. Full line: classical 
model; dashed line: both relativistic models. 
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Fig. 6 Eccentricity evolution of the outer planetary system. Full line: classical 
model C; dashed line: relativistic model Rs ; dashed dotted line: relativistic model 

Rail- 
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Fig. 7 Inclination evolution of the outer planetary system. Full line: classical model 
C; dashed line: relativistic model Rs ; dashed dotted line: relativistic model R ;;. 
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semimajor axis (AU) 
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semimajor axis (AU) 



Fig. 8 Maximum differences between cfassic and reiativistic secuiar dynamics gen- 
erated by the giant pianets as function of the initial (a, q) values in a numerical 
integration for 10 5 years. All test particles have initial i = 20° , S2 = 0°, u = 60° . 
The extreme changes observed in (e, i) at a ~ 0.4 are due to a strong Kozai mech- 
anism (see Fig. 9). These differences increase when considering the secular effects 
of the entire planetary system. 
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Fig. 9 Time evolution of a fictitious particle inside the secular Kozai mechanism 
generated by the giant planets. Full line: classical model C; dashed line: rela- 
tivistic model Rs. In spite of considerable variations in e and i the parameter 
H — vl — e 2 cosi remains nearly constant in both models. Relativistic terms di- 
minishes the amplitude of the oscillations in (e, i). 
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Fig. 10 Asteroid 2004 XY60. At the beginning the asteroid is in the 6:5 resonance 
with Venus. After 13000 years the strength of the resonance drops and both the 
classic and relativistic particles are captured in different high order mean motion 
resonances with Earth. 
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Fig. 11 Asteroid 2000 LK. An example where relativistic effects plus chaotic dy- 
namics generate very different evolutions. Full line: classical model C; dashed line: 
relativistic model Rg. The relativistic simulation ends up colliding with the Sun. 
In this case the parameter H is not conserved. 
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Fig. 12 Asteroid 2003 CP20. An example of evolution without encounters with 
the planets. Full line: classical model C; dashed line: relativistic model R5. The 
asteroid is captured in resonance 27:28 with Venus and follows the effects of the 
Kozai mechanism. 
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Asteroid 2003 CP20. The critical angle of the resonant motion is a 
+ 28A — vj according to the classic and relativistic models. 



